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La 2 03 Fe 2 Se 2 , a Mott insulator on the brink of orbital-selective metalization 

Gianluca Giovannetti, 1-2 Luca de’ Medici, 3 Markus Aichhorn, 4 and Massimo Capone 2-1 

1 CNR-IOM-Democritos National Simulation Centre, 

UoS Trieste-SISSA, Via Bonomea 265, 1-34136, Trieste, Italy 
2 International School for Advanced Studies (SISSA), Via Bonomea 265, 1-34136, Trieste, Italy 
3 European Synchrotron Radiation Facility, BP 220, F-38043 Grenoble Cedex 9, France 
4 Institute of Theoretical and Computational Physics, TU Graz, Petersgasse 16, Graz, Austria 

We show that the insulating character of the iron-selenide La 203 Fe 2 Se 2 can be explained in 
terms of Mott localization in sharp contrast with the metallic behavior of FeSe and other parent 
parent compounds of iron superconductors. We demonstrate that the key ingredient that makes 
La 2 C> 3 Fe 2 Se 2 a Mott insulator, rather than a correlated metal dominated by the Hund’s coupling is 
the enhanced crystal-field splitting, accompanied by a smaller orbital-resolved kinetic energy. The 
strong deviation from orbital degeneracy introduced by the crystal-field splitting also pushes this 
materials close to an orbital-selective Mott transition. We predict that either doping or uniaxial 
external pressure can drive the material into an orbital-selective Mott state, where only one or few 
orbitals are metallized while the others remain insulating. 

PACS numbers: 71.30.+h, 71.10.Fd, 71.27.+a 


I. INTRODUCTION 

The link between high-temperature superconductiv¬ 
ity and strong electron-electron correlations has been 
forged and strengthened by decades of investigation in 
the copper-based superconductors (cuprates). In this 
light, the debate about the strength and the role of elec¬ 
tron correlations in iron-based superconductors (FeSC) 
maintains a crucial importance. The overall phenomenol¬ 
ogy of these materials does not provide a self-evident an¬ 
swer. In these materials superconductivity appears dop¬ 
ing a metallic spin-density-wave parent compound, rather 
than the Mott insulator of the cuprates, and, while the 
metallic state is highly incoherent, the standard finger¬ 
prints of strong correlations, like the Hubbard bands, are 
not universally observed. 

Furthermore non-perturbative studies of the effect of 
the interactions, mainly based on Dynamical Mean-Field 
Theory (DMFT) and related methods, have highlighted 
a novel behavior, in which the Hund’s coupling Jy, to play 
a major role in determining the degree of correlations^. 
The electrons in the d-orbitals (the parent compounds 
have a nominal filling of six electrons for each iron atom) 
are strongly correlated, as measured by the bad-metallic 
behavior with small coherence scales shown in many ex¬ 
periments, but the Hubbard repulsion U is substantially 
smaller than the critical value for the Mott transition^. 
This regime is often labeled as a “Hund’s metal” and 
displays characteristic properties^ including a remark¬ 
able tendency to “orbital selectivity”, i.e., to a neat dif¬ 
ferentiation in the degree of correlation of the different 
orbitals, leading to the simultaneous presence of weakly 
and strongly correlated electrons^! - — . In Ref— it has 
been shown that the Hund’s coupling decouples the or¬ 
bitals quenching the inter-orbital fluctuations and that 
this leads to a picture of five doped single-band Mott in¬ 
sulators. Consequently the degree of correlation is con¬ 
trolled by the distance of each individual orbital from the 


half-filled d 5 configuration. The relevance of the d 5 Mott 
phase in the phase diagram of models for iron supercon¬ 
ductors has been observed also in 

While no iron-based material can be doped with one 
hole per iron site, isostructural materials where iron is re¬ 
placed with manganese are characterized by the d 5 half- 
filled configuration and they are indeed antiferromagnetic 
Mott insulators— r—. However, the evidence for strong- 
correlation physics in the FeSC does can not rely on the 
presence of actual Mott states directly connected with 
the superconducting compounds. Indeed no Mott insu¬ 
lator exists in the 122 family (BaFe 2 As 2 or isoelectronic 
compounds doped either with holes and electrons), the 
1111 family which originates doping LaFeAsO, the 11 
selenides FeSe and FeTe, and LiFeAS. While this experi¬ 
mental fact is indeed completely compatible with the sce¬ 
nario based on the Hund’s coupling, it may cast doubts 
on the whole relevance of correlations. 

A close relative of FeSC with insulating behavior is 
La 2 03 Fe 2 Se 2 . This system is based on a square lattice of 
Fe ions with nominal valence 2+ as in all of iron pnictides 
and chalcogenides. However, the resistivity as a function 
of temperature shows clearly an insulating behavior with 
an activation energy gap of about 0.19eV— On the ba¬ 
sis of this experimental evidence combined with Density- 
Functional Theory (DFT) supplemented by mean-field 
treatment of the Hubbard U (DFT+U) in the ordered 
magnetic phase it has been argued that La 2 03 Fe 2 Se 2 
is in a Mott state with low-temperature antiferromag¬ 
netic ordering.*! 5 This evidence raises the natural ques¬ 
tion about the reason why the Hund’s metallic behav¬ 
ior of undoped FeSC is replaced with a Mott state in 
this compound. A simple bandwidth reduction is a very 
unlikely answer, since one of the defining properties of 
the Hund’s metal regime is that the critical value for a 
Mott transition is pushed to very large values, signifi¬ 
cantly far from the experimental estimates.— This calls 
for the identification of specific aspects of the bandstruc- 
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ture which make the present material insulating. 

In this work we investigate the differences of the elec¬ 
tronic structure of La 2 C> 3 Fe 2 Se 2 with respect to parent 
compounds of the FeSC using the reference case of FeSe. 
We identify in a strong lifting of the orbital degeneracy 
the main difference of the material under consideration 
with respect to most parent compounds of iron-based su¬ 
perconductors. This leads to a marked tendency towards 
orbital-selective Mott transition and opens the way to 
a complete Mott localization when combined with the 
reduced kinetic energy of La 20 3 Fe 2 Se 2 with respect to 
FeSe. 

The paper is organized as follows. In Sec. II we present 
the electronic structure which emerges from Density- 
Functional Theory calculations, while In Sec. Ill we 
consider the inclusion of electron-electron interactions. 
The section is organized in two subsections dedicated 
respectively to the results obtained within the slave- 
boson mean-field scheme and the more accurate Dynami¬ 
cal Mean-Field Theory. Sec. IV contains our conclusions. 


II. DENSITY-FUNCTIONAL THEORY 
BANDSTRUCTURE 

We start our investigation determining the bandstruc- 
ture of La 2 C> 3 Fe 2 Se 2 by means of DFT in the frame¬ 
work of the generalized gradient approximation (using 
the Perdew-Burke-Ernzerhof (PBE) functional^) for the 
tetragonal unit cell of La20 3 Fe 2 Se2^ using Quantum 
Espresso^ and Wien2Ki£. In Fig. Q]we show the density 
of states projected on the different atoms. 

PBE calculations clearly lead to a metallic solution 
with a sizable spectral weight at the Fermi level. This 
low-energy contribution to the spectral density is domi¬ 
nated by the bands arising from Fe 3d electrons, which 
are very weakly entangled with oxygen and lanthanum 
bands lying in energy windows far from the Fermi level. 
The 3d bands have an overall width of 3.2 eV, consid¬ 
erably narrower than the 4.6 eV of the same bands in 
FeSoil 

In order to include the on-site Coulomb interaction pa¬ 
rameterized by the Hubbard U and the Hund’s coupling 
Jh we compute maximally localized Wannier orbitals^S 
for the pure 3d Fe orbitals built from the iron bands in 
the energy range between -2 and 1.2 eV. The properties 
of these orbitals will also provide us with important in¬ 
formation that will help us to rationalize the behavior of 
this materials. 

The on-site energies of the Wannier orbitals reflect 
indeed an important difference between the material 
under consideration and the parent compounds of the 
iron-based superconductors. In La 2 C> 3 Fe 2 Se 2 , the Fe 
ions are surrounded by two nearby oxygen anions and 
four distant selenium anions. The additional oxygen 
ions result in a lower symmetry in La 20 3 Fe 2 Se 2 as 
compared to FeSe and other FeSC. This means that 
the local problem is no longer diagonal in the stan¬ 



-8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 

E(eV) 


FIG. 1. (Color online) Density of states of La20sFe2Se2 cal¬ 
culated within PBE and projected on the different atomic 
species. The zero energy is set at the Fermi level. 


dard cubic basis defined by the e g and t 2 g orbitals. 
We can obviously diagonalize the local Hamiltonian for 
La20 3 Fe2Se2 by means of a unitary transformation. The 
resulting orbitals, which arc linear combinations of the 
(3z2 — r2,xz,yz,x2 — y2,xy ) orbitals^i that we label 
as (1,2,3,4,5) have in our calculations on-site energy 
of (-0.660, -0.627, -0.397,0.183,0.543) eV, respectively. 
This leads to a total crystal-field splitting of 1.1 eV, much 
larger than the value for FeSe, which we estimate in 
0.48 eV within analogous PBE calculations. 

Comparing the orbital-resolved density of states for 
La 2 0 3 Fe2Se 2 and FeSe (see Fig. EH we can visualize the 
different crystal-field splitting and emphasize another re¬ 
lated effect: in FeSe the different orbitals have a similar 
width and lie more or less in the same energy range. In 
contrast La 20 3 Fe 2 Se 2 features a set of narrow orbitals 
significantly shifted in energy relative to each other. That 
means that the effective band-width reduction for each 
individual orbital is even larger than the factor estimated 
from the total band width. Interestingly, even in FeSe 
and other standard FeSC it has been showed that, ar¬ 
bitrarily neglecting interorbital hybridizations turns the 
broad orbitals into more localized objects, similar to what 
we find for La20 3 Fe2Se2. Therefore we can link the pecu¬ 
liar orbital-resolved density of states of the latter material 
with the crystal-field splitting that reduces the effect of 
the interorbital hybridization. 

In the following of the paper we investigate how the dif¬ 
ference in the single-particle spectra reflects on the effects 
of electron-electron correlations, keeping in mind the role 
of the Hund’s coupling that we anticipated above. The 
results picture can be quite rich because of the non-trivial 
interplay between the “Hund’s physics” and the two main 
differences we highlighted, a reduced kinetic energy and 
an increased crystal-field splitting. In particular, while a 
reduced kinetic energy simply leads to effectively larger 
Coulomb terms, the crystal-field splitting competes with 
the Hund’s coupling, as the latter tries to spread the elec- 
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trons among the different orbitals to maximize the total 
spin, while a large crystal-field splitting obviously favors 
an unbalanced population with large occupation of the 
low-lying orbitals. 


III. EFFECT OF ELECTRON-ELECTRON 
CORRELATIONS 


To understand the role of electron-electron interactions 
in turning La 2 C> 3 Fe 2 Se 2 insulating we consider two dif¬ 
ferent approaches to treat the short-range interactions. 
We start from a Slave-Spins mean-field (SSMF) theory^, 
which allows for a computationally inexpensive and fast 
survey of the phase diagram and it is expected to cap¬ 
ture the main physics as long as the system remains in 
a Fermi-liquid stato^. Then we move to the more ac¬ 
curate DFT+Dynamical Mean-Field Theory (DMFT)— 
method, which treats exactly the local quantum dynam¬ 
ics mapping the lattice model onto an impurity embed¬ 
ded in a self-consistent bath. As impurity solver we em¬ 
ploy mainly Exact Diagonalization (ED) j? - i — We ver¬ 
ified for selected parameter the excellent agreement of 
ED with the continuous-time Quantum Monte Carlo so¬ 
lution of the impurity model implemented in the TRIQS 
toolkit^ - — . 

For the interacting part of the Hamiltonian we use a 
Kanamori form parameterized by U and Jh according to 


Hint — U ^ ^ ’hlima'^iTrur' T U ^ ^ 

i,m i,m,m 

+ u" E W'ima'ft'im' a 


nmff'Hm'u' 


+ dtm-tdtmldim'tdim'l} 


(1) 


where dq mcr is the destruction operator of an electron 
of spin a at site i in orbital m, and ni ma = dJ^dima, 
U and U' = U — 2 Jh, U = U — 3Jh are intra- and 
inter-orbital repulsions and Jh is the Hund’s coupling. 
In the absence of estimates of U and Jh for this mate¬ 
rial, we use the constrained random-phase-approximation 
(cRPA) for FeSe (U = 4.2 eV and Jh = 0.504eV) which 
are not expected to differ in a critical way* 22 Note that 
these interaction values are given in Hubbard-Kanamori 
notation. The use of the same parameters also allows 
us to highlight the role of the material-specific properties 
(kinetic energy and crystal-field splitting) in determining 
the low-energy properties of the system. In all the calcu¬ 
lations we consider paramagnetic solution that preserve 
also the orbital symmetry. 
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FIG. 2. (Color online) Orbital-resolved density of states for 
La 2 0sFe2Se2 (bottom panel) and FeSe (top panel). The zero 
energy is set at the Fermi level. 


A. Slave-spin Mean-Field Approximation 

In Fig. [3] we show the DFT+SSMF orbital resolved 
densities n m and quasiparticle weights Z m as functions of 
Coulomb parameter U at fixed ratio Jh/U= 0.2. Within 
this approximation spin-flip and pair-hopping terms (last 
two lines of Eq. Jl])) are not included. Z m is the measure 
of the metallic character of each orbital. A vanishing of 
Z m for a given orbital implies that the carriers with that 
character are reaching a Mott localization. A Mott insu¬ 
lating state is reached when all the Z m 's vanish, while a 
situation with coexisting finite and zero Z n J s would cor¬ 
respond to an orbital-selective Mott transition (OSMT). 

For small interactions, up to U ~ 1.5eH, the orbital 
densities range from 0.22 to 0.8 and the Z m are close 
to one. Increasing U, the Z m decrease and depart from 
each other while the orbital populations deviate substan¬ 
tially from the non-interacting values. More precisely, 
the orbitals labeled as 4, 3 and 5 move, one after the 
other, towards a half-filled configuration (ro m = 0.5) and 
they become insulating ( Z m = 0) at different interaction 
strengths. In other words, the system shows a series of 
OSMT’s in which the different orbitals become localized 
independently on the behavior of the others. The two 
remaining orbitals (1 and 2) instead remain metallic for 
a larger range of interactions and they simultaneously 
become insulating for U c = 3.6eH, smaller than the es¬ 
timated value for FeSe. At this critical interaction also 
orbital 2 becomes half-filled, while orbital 1 is completely 
filled. The “order” in which the different orbitals undergo 
a Mott transition is a consequence of the individual band¬ 
width, which is smaller for orbitals 4 and 3 and of the 
initial orbital population. 

The series of transitions that we described is in stark 
contrast with the results for FeSe^ and for the other par¬ 
ent compounds^^, where no OSMT occurs despite a 
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marked differentiation between the Z m of the different or¬ 
bitals. In this case a full Mott transition in which all the 
orbitals become simultaneously localized takes place for 
values of the interactions which are much larger than rea¬ 
sonable estimates for this material. In particular, for the 
cRPA values (for FeSe) introduced above, La 2 C> 3 Fe 2 Se 2 
is a Mott insulator, while FeSe is a Hund’s metal with 
small and orbital-differentiated quasiparticle weights. 

It is useful to notice that, increasing the value of Jh the 
critical coupling for an OSMT is reduced, while the crit¬ 
ical U for a full Mott transition increases^. This result, 
together with the whole picture we have drawn, is in per¬ 
fect agreement with model calculations in which orbitals 
with the same bandwidth have been shifted leading to an 
OSMT— 

Interestingly, while the full Mott localization requires 
a commensurate filling, the OSMT that we observed sur¬ 
vive also doping the d 6 Mott insulator. In the inset of 
Fig. [3] we show the evolution of the quasiparticle weight 
as a function of U for a d 5 " 5 configuration, where we 
doped one hole every two iron sites. Comparing with the 
d 6 case, we see that orbitals 3, 4 and 5 still undergo the 
series of OSMT’s we described above, and in particular 
the three orbitals are insulating for the U values repre¬ 
sentative of La 2 C> 3 Fe 2 Se 2 , while orbitals 1 and 2 become 
metallic after doping. This can be explained in terms of 
a schematic general picture of OSMT—. For sizable J/U 
the global half-filled configuration is in general a Mott in¬ 
sulator at quite smaller U than for the d 6 case. In this d 5 
Mott insulating phase each orbital opens an independent 
gap, thus being half-filled when the chemical potential 
falls within its gap. The gaps can have different widths 
and position depending on the crystal-field splitting and 
the bandwidth of their orbital. Doping of a Mott insula¬ 
tor occurs essentially when the chemical potential moves 
out of the gap, and thus in this situation this happens 
for each orbital at a separate value of the chemical poten¬ 
tial. The selective doping of one or more orbitals leads 
therefore to an OSMT. Our results show that the band- 
structure of La 203 Fe 2 Se 2 is such that when reaching a 
total filling of 6 electrons only two orbitals (1,2) have ac¬ 
tually been doped compared to the half-filled case while 
the rest remain half-filled, have an open gap, and are 
thus insulating. This does not depend on the filling be¬ 
ing exactly 6 or slightly above or below, thus doping the 
d 6 configuration tunes indeed the Mott transition for or¬ 
bitals 1 and 2 but does not alter the selective insulating 
behavior for orbitals 3,4,5. 

Our results strongly indicate that the main anomaly 
of La 2 0 sFe 2 Se 2 with respect to the parent compounds of 
FeSC is the strong enhancement of the crystal-field split¬ 
ting with respect to iron-based superconductors, which 
conjures with the Hund’s coupling to decouple the dif¬ 
ferent orbitals, leading to an orbital selective localization 
rather than to a global Mott transition for the six elec¬ 
trons. The overall reduction of the kinetic energy with 
respect to FeSe plays a more quantitative role, reducing 
the actual critical value of U. 
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FIG. 3. (Color online) Orbitally resolved density (n m ) and 
quasiparticle weights ( Z m ) calculated within DFT + SSMF 
for La20aFe2Se2 at ratio Jh/U = 0.20. The inset shows Z m 
for a filling of d 5 ' 5 electrons per iron (one hole every two iron 
sites). 


B. Dynamical Mean-Field Theory 

Our DFT+SSMF treatment is however limited to 
metallic solutions and it can not help us to characterize 
the Mott insulating solution for U > U c . A full descrip¬ 
tion of both low- and high-energy features can instead 
be obtained by means of DFT+DMFT. As mentioned 
above, we mainly use and ED solution for the impurity 
model, which also allows us to include pair-hopping and 
spin-flip terms, as opposed to DFT+SSMF. 

We consider N a = 15 orbitals in total, with 5 impurity 
orbitals and 2 bath degrees of freedom connected to each 
impurity orbital. The storage requirements for the ma¬ 
trix hamiltonian and Lanczos vector, since the Hilbert 
space has dimension 41 409 225, are solved by splitting 
them over parallel processors. The impurity solver is di¬ 
agonalized by a parallel Arnoldi algorithm^ using the 
symmetry with respect to the inversion of up and down 
particles. 

In Fig. [3] we show the dynamical information obtained 
for the cRPA values U = 4.2 eV and Jh = 0.504 eV. 
We present orbital-resolved Green’s functions on the 
imaginary-frequency axis (a)), which is the most direct 
product of the calculation, and the real-frequency spec¬ 
tral function (b)), which provides a more obvious physical 
content. The imaginary part of the Matsubara Green’s 
functions smoothly extrapolate to zero in the limit of 
small frequency. This clearly confirms that the material 
is insulating also within the more accurate DFT+DMFT 
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method. The imaginary part of the self-energy diverges 
only for orbital 5, for which the chemical potential hap¬ 
pens to lie in the middle of the gap (See also the spectral 
function in panel b)). For the other orbitals the chemical 
potential lies far from the center of the gap and, conse¬ 
quently, the imaginary part of the self-energy does not 
diverge. Finally, a comparison of the orbital-projected 
spectral density with the DFT results of Fig. [Hallows us 
to visualize the formation of high-energy spectral weight 
due to strong correlations and to estimate the Mott gap 
A Mott — 1 eV . The Mott insulating state of La20 3 Fe 2 Se2 
is characterized by a high-spin configuration with S = 2 
as a consequence of the Hund’s coupling. The orbital 
populations are in agreement with the DFT+SSMF re¬ 
sults, confirming the picture of four half-filled orbitals 
and one full orbital. Therefore the whole scenario ob¬ 
tained within DFT+SSMF survives when we move to the 
full DMFT treatment of correlations, suggesting that also 
the sequence of OSMT transition preceding the full Mott 
localization is a real feature of the present material. The 
disagreement between the theoretical value of the gap 
and experiments may be attributed to an overestimate of 
U and Jh , or even in their ratio. However, as mentioned 
above, the choice to use the same values of FeSe allowed 
us to highlight the difference between La 203 Fe 2 Se 2 and 
the parent compounds of standard FeSC. 

All the calculations we presented are limited to a 
paramagnetic solution without broken symmetry in any 
channel, demonstrating the intrinsic Mott character of 
the insulating state. However, at low temperature this 
state may be unstable towards different ordering as mag¬ 
netic and orbitals one s 33 ’ 34 as it is indeed suggested by 
DFT+U calculations^. 

To confirm the reliability of our results, we com¬ 
pared with a different DFT+DMFT implementation 
based on the TRIQS package^ 7- — In these calculation 
we use projective Wannier function techniques, as well 
as continuous-time quantum Monte Carlo as impurity 
solver. In addition, we considered also the Se p states 
for the construction of the Wannier functions, in order 
to check the stability of the insulating state. For the in¬ 
teractions we use also cRPA values for FeSe, and a Slater 
expansion of the interaction hamiltonian*. 

The results confirms the precedent picture. The sys¬ 
tem is in a Mott state, with four orbitals close to half- 
filling, and one completely filled. The insulating state 
appears therefore extremely solid, as it does not depend 
on the details of the computational scheme that is ap¬ 
plied. 


IV. CONCLUSIONS 

In conclusion, we use different approaches combin¬ 
ing DFT with strong correlation physics to study the 
electronic properties of the new synthesized compound 
La 2 03 Fe 2 Se 2 , an insulating material that shares the same 
electron count with the metallic parent compounds of 
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FIG. 4. (Color online) a) Orbital-resolved imaginary part of 
local Green’s function (G) as a function of the Matsubara 
imaginary frequency, b) spectral density (A) as a function 
of the real frequency computed within DFT + DMFT for 
La203Fe2Se2. 


iron superconductors. Our calculations clearly demon¬ 
strate that the material is a Mott insulator. A detailed 
understanding of the electronic structure of La 203 Fe 2 Se 2 
compared to the one of FeSe highlights the difference of 
this class of compounds with ordinary Fe-based super¬ 
conducting materials. Using the same interaction values 
for the two materials, we demonstrate that the insulating 
behavior of La 2 03 Fe 2 Se 2 is not simply due to a larger 
interaction strength, but it is related to a more funda¬ 
mental difference in the electronic structure. 

We identify the main reason for the insulating behav¬ 
ior in the enhancement of the crystal-field splitting due 
to the position of the oxygen atoms. The larger crystal- 
field leads to a reduced overlap between the density of 
states with different orbital character, which triggers a 
series of successive Mott transitions in which the dif¬ 
ferent orbitals become insulating one after the other as 
the interaction grows. The insulator is indeed character¬ 
ized by four singly occupied orbitals and a fully occupied 
one. This condition is quite different from iron-based 
superconductors, in which a marked orbital differentia¬ 
tion does not however lead to individual Mott insulating 
behavior. 

As a consequence, we predict that La20 3 Fe 2 Se2 should 
give rise to a series of OSMT under pressure and we ex¬ 
pect that a similar behavior should take place as a func¬ 
tion of doping. A slightly doped material is expected 
to show the coexistence between localized orbitals and 
delocalized orbitals hosting the extra charges. 

Note Added. During completion of this work, we be¬ 
came aware of a related work by B. Freelon et Al., which 
is complementary to ours and gives comparable results. 
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